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^ Abstract 
CO' 

The transport properties of Brownian ratchet was studied in the presence of stochastic intensity 
noise (SIN) in both overdamped and underdamped regimes. In the overdamped case, analytical 
solution using the matrix continued fraction method revealed the existence of a maximum current 
when the noise intensity fluctuates on intermediate time scale regions. Similar effects were observed 
for the underdamped case by Monte Carlo simulations. The optimal time-correlation for the 
Brownian transport coincided with the experimentally observed time-correlation of the extrinsic 
noise in Esherichia coli gene expression and implied the importance of environmental noise for 
molecular mechanisms. 
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1. INTRODUCTION 



Noise-induced phenomena are attracting much attention not only in engineering but 
also in molecular biology. Counter-intuitively, noise can enhance system performance by 
increasing transmission and synchronization of information through stochastic resonance 



6J and noise- induced synchronization |7H10|. 
For example, molecular systems function efficiently in nano-scale environments under 
multi-scale noise through thermal and other environmental fluctuations llNl4|. This ef- 



ficiency should not be interpreted by assuming steady-state or Gaussian distribution. A 
recent single-cell observation of Escherichia coli revealed that the protein copy number does 



not obey the gamma distribution 



15], and its stationary distribution can be approximated 



by superstatistics (i.e., the superposition of multiple statistical models; see Equation [T] and 
Discussion). Chabot et al. revealed that the cellular variability in gene expressions comes 



from temporal (periodic) noise which is related to circadian oscillation 



In accordance 



with such experimental observations, theoretical studies also conclude that biochemical noise 
is not Gaussian to facilitate enhanced functionality [xs| - Both experimental and theoret- 
ical approaches suggest fundamental roles of noise-enhanced phenomena to render efficient 
molecular systems. 



In this paper, we investigate the efficiency of Brownian motor |19l . [20| (or a ratchet 
transport) under a noisy, nonequilibrium state. It is known that the violation of detailed 
balance induces a transport effect, which is a model for many biological mechanisms including 



ion pumps 



21 



23] that use ATP for transport (In theory, however, the conformational 



fluctuation of such pumps can facilitate transport even without ATP 



Brownian transport has been intensively studied including in mass separation |25J, elec- 
tron transport in a q „a„ tum ratchet fl — of ato mS !„ optical traps 0, a "ndo m 
walker model 28L and non-Gaussian noise models 29|, |3Gj (for more details, see compre- 
hensive reviews 20, 3j|). One of the most popular ratchet models is the correlation ratchet 
in which Brownian particles in a ratchet potential are driven by the addition of white and 
colored noise. The model studied here is a ratchet driven by the multiplication of white 
and colored noise. Let us introduce our model more formally. 

Superstatistics with tem por al and/or spatial fluctuations 32j-|35] is used to explain non- 



Gaussian distributions 



36 



37] in applied physics. This concept is also seen in stochastic 
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processes in which noise fluctuation is treated in a static way [33l. |38M43| . The superstatistical 
stochastic process calculates a stationary distribution P s t(x) by assuming that the noise 
fluctuates over a long time scale (i.e., very slowly): 



Pst(x) 



dDP 8t {x\D)P{D), 



(1) 



where P s t(x\D) is the stationary distribution given the noise intensity D and P{D) is the 
distribution of D. Equation [I] fits in the Bayesian framework by considering P s t(x) and 
P{D) as posterior and prior distributions, respectively. 

We model temporal noise-intensity fluctuation dynamically and modulated the inten- 
sity of white Gaussian noise by the Ornstein-Uhlenbeck process in overdamped Langevin 



equations 



44 



45|: 



x = -V'(x) + sUt), 

a = -j(s -a) + yfy£ a (t). 



(2) 
(3) 



Here, V(x) is the potential, a is the mean of the Ornstein-Uhlenbeck process, 7 is the 
relaxation rate (7 > 0), and £ x (t) and £ s (t) are white Gaussian noise with the correlation: 



(UW)) = 2DJ(t-lf), 
(&(*)6(0> = 2D s 5(t-t'). 



(4) 
(5) 



We call the term s(t)£ x (t), the stochastic intensity noise (SIN) because the noise intensity is 
governed by a stochastic process. SIN is the multiplicative term of white and colored noise, 
and qualitatively different from white noise: it is in nonequilibrium. 

In the context of Brownian transport, Reimann et al. 46| first studied the transport effect 
with sinusoidal noise-intensity modulation. Our work differs from this and succeeding studies 



46^481] 



m our 



that employed a discrete dichotomous noise or a deterministic periodic signal 
model, fluctuations are governed by a continuous stochastic process (the Ornstein-Uhlenbeck 
process). There exist similar models. Borromeo et al. studied a current generated by two 
symmetric colored noises, the Ornstein-Uhlenbeck noise and its time-delayed version, and 
observed the Maxwell's daemon-like phenomenon 49(. Morgado et al. investigated temporal 
heterogeneity in Poisson mechanism [5^] . Our model focuses on a multiplicative rather than 
an additive effect, because biological phenomena are governed by multiplication (also see 
Discussion). This distinction highlights the importance of SIN-induced transport. 



In our calculations, we investigated the effect of four controllable parameters on the 
current: 

• 7 (the relaxation rate in Equation [7]) , 

• Q (the effective noise intensity in Equation [11]) , 

• p (the squared variation coefficient in Equation [T5J) and 

• /i (the scaled mass in Equation [H]). 

Although D x (the noise-intensity in Equation [6]) is also a controllable parameter, we kept 
it constant (D x = 1) throughout the paper. The squared variation coefficient p, which is 
generally defined as the ratio of the squared mean to the variance, characterizes the deviation 
of the Gaussian noise distribution by the kurtosis. The overdamped case is calculated 
using the matrix continued fraction method (MCFM) and Monte Carlo (MC) simulations. 
The calculations reveal that the current is a maximum at adequate 7 and Q. This result 



concurs with resonant activation (RA) [5l| and noise-enhanced stability (NES) [52|, |53[ in 
the escape problem. Our main result is the enhanced transport capability in intermediate 
time-correlation regions. This has an important biological implication. The time-correlation 
of extrinsic noise (i.e., environmental fluctuation) with gene expression in E. coli is on the 



order of cell cycle length |54|. By fitting the time-scale of our model to the E. coli model in 
Ref. {55], we found that the time-correlation of extrinsic noise coincides with the region of 
enhanced transport. It is a theoretical backup of the tactful exploitation of environmental 
fluctuation by biological organisms. 

2. METHODS 

2.1. Brownian Particles 

Brownian particles are subject to noise-intensity fluctuations represented by 

Hx = -x-V'(x)+s£ x (t), (6) 
s = _ 7 ( s _ Q! ) + ^ a (t) j (7) 

where the scaled mass p is introduced into Equation [2] for a study of a mass separation effect, 
and meanings of V(x) are the same as in Equation [2j We used the same ratchet potential 
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FIG. 1: (Online version in color) A ratchet potential of Equation [HI The dashed and solid arrows 
indicate the directions of the current in the correlation ratchet (driven by additive white and colored 
noise) and the ratchet (driven by SIN), respectively. 

(b) p = 0.1 



(a) p = 0.01 



(c) P =l 



= 100 




FIG. 2: (Online version in color) (a)-(d) Time courses of SIN generated by MC simulations for 
four values of p (squared variation coefficient given by Equation fT2|) : p = (a) 0.01, (b) 0.1, (c) 1, 
and (d) 100. (e)-(h) Histograms of SIN for the four values of p and their sample kurtosis k. In 
(a)-(h), we varied p while keeping the effective intensity Q = D X (D S + a 2 ) constant. 



function as used in previous studies [56|, [57] 

V(x) ' 



sin(27rx) + - sin(47rx) > + Fx, 

Z7T 4 



(8) 



with a periodicity V(x+1) = V{x) (F = 0), where F is the load. Figured] shows the potential 
with no load (F = 0), where the dashed and solid arrows indicate the normal current 
direction for the correlation ratchet [sol, 57 1 and for the ratchet driven by SIN, respectively. 
The current direction for the SIN case is identical to that in Ref. ^7| in which the noise 
intensity is modulated by a random dichotomous process. In Equation [3, the relaxation rate 
7 denotes the inverse of a time scale (time-correlation). When the noise intensity fluctuates 
with a longer time scale (7 — > 0), systems driven by SIN locally equilibrates and hence a 
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current is not generated in accordance with the second law of thermodynamics. Likewise, 
SIN reduces to white noise (with a noise intensity Q) when the noise-intensity fluctuates 



very rapidly (7 — > 00) [44|, which also indicates that the current vanishes. 

For Equation the stationary distribution P s t(s) of the intensity-modulation term s is 
given by 

PJs) = —=t= exp I l^(s-a) 2 \. (9) 



A calculation of the correlation function of SIN yields [45 1 

(s(t)£ x (tW)Ut')) = 2Q6(t - 0, (10) 

with 

g = J D x ( J D s + a 2 ), (11) 
where Q expresses the effective noise intensity and D x and D s are the noise intensities of 



^x(t) and £ s (i), respectively [Equations 0] and 
coefficient p of the noise-intensity fluctuations 



. Here, we introduce the squared variation 



P = D s /a 2 , (12) 

which denotes the squared ratio of the standard deviation to the mean of Equation [9] 
Figures Eta)-(d) show trajectories of SIN with p = 0.01, 0.1, 1, and 100, respectively. All 
share the same effective noise intensity Q (the same variance). For p = 0, SIN reduces to 
white Gaussian noise with a noise intensity Q = D x a 2 . On increasing p (Figures [2](e)-(h)), 
the noise-intensity fluctuations become larger and the distribution of SIN deviates from the 
Gaussian distribution. To quantify the deviation, we calculate the kurtosis of SIN, which 
is a measure of heavy tails in probability density functions. The kurtosis k of SIN (i.e., 
s(t)^ x (t)) is given by (Appendix [AT) 

{{s(t)Ut)V) n 6 



K 



9- n , N2 (0<p<oo), (13) 



[{ S (t)ut)} 2 ) 2 ( i+ py 

where k depends only on the squared variation coefficient p having a crucial effect on the 
statistical properties of SIN. Equation [13] shows that the kurtosis is 3 for p = and k 
increases with increasing p, giving a flatter distribution. Equation [13] is plotted by the solid 
curve in Figure [3] in which the filled circles denote the kurtosis calculated by MC simulations. 
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FIG. 3: (Online version in color) Kurtosis k of SIN as a function of p. The solid line indicates 
Equation [13] and the filled circles indicate the sample kurtosis in Figures [2{e)-(h). 

3. RESULTS 

In Equations |6] and our model has four parameters (7, a, D x , and D s ) in terms of noise 
properties. Instead of these parameters, we adopt 7, p (the squared variation coefficient 
in Equation [T2l) . Q (the effective noise intensity in Equation [TTi) . and D x as controllable 
parameters in model calculations, where Q and p are defined by Equations [TT] and [£2l 
respectively. The new parameters (7, p, Q, and D x ) specify the system identically since a 
and D s are uniquely determined by D x , Q, and p: 



Performing model calculations based on the MCFM and MC simulations, we studied the 
dependence of the current on 7, p, and Q in the overdamped case (D x is set to D x — 1 
throughout the paper). In the underdamped case, we additionally investigated the depen- 
dence of the current on the scaled mass p (Equation [6]) with MC simulations. MC simulations 
were performed with the Euler-forward method with a time resolution of At = 10 -4 (for 
details of the method, see Ref. [581]) . 

3.1. Overdamped case 




(14) 



We first calculated the current J for the overdamped case (p = in Equation [6]). For 
t — > 00, the stationary distribution P st (x, s) of (x, s) has to satisfy the stationary Fokker- 

7 



Planck equation (FPE) 

h FP P st (x,s) = 0, (15) 

where Lpp is an FPE operator: 

Lfp = |y. (l) + Dm <jLf +1 {| (> -<•)+!>.£}. (16) 

Since Equation [15] is not written in terms of potential forms, we cannot calculate the station- 
ary distribution in a closed form. Consequently, we used the MCFM to solve Equation [T5] 
which expands P st (x,s) in terms of a complete orthonormal set. The MCFM is a common 



technique for stochastic processes and it is widely used to solve FPEs (see Ref. 58] and 
references therein). Considering the periodicity of the potential [V'(x + 1) = V(x)] and 
domain [x G (— oo, oo) and s G (— oo, oo)], we expanded the stationary distribution P s t(x, s) 
in a Fourier series for x and the Hermite function for s: 

Pst(x, s) = ip (s) 22 ^2C kjn exp(27ckix)^ n ( y s). (17) 

k=-M k n=0 

Here, Ck, n are expansion coefficients, Mk and M n are truncation numbers on which the 
precision of obtained solutions depends, and if) n (s) is the Hermite function satisfying the 
orthonormality relation (^) n i (s)t/} n (s)) = 5 n i n \ 

1/4 



(ak) vS ff " ( "' CTp (-^) " = &" a) ' (18) 



where H n (s) is the nth Hermite polynomial. Multiplying Equation [15] by 
exp(— 27ik'ix)ipn'(s)/ip ( y s) and integrating over x and s, we obtain a linear algebraic equa- 
tion in terms of Cfc jn , which can be solved by the MCFM (Appendix [B]). The current J is 
calculated using 

J = (v) = J dx(^j dsJ x (x,s)j, (19) 

where J x (x,s) is the probability current in the x direction due to the continuity equation 
d t P + d x J x + d s J s = 0: 

J*(x, s) = (-V'{x) - D x s 2 ^j P st {x, s). (20) 
By substituting Equation [17] into Equation [19] the current J can be expressed in terms of 

Ck,n- 

J = ~{C-i, + C lfi ) - \ (C_ 2 ,o + C 2i0 ) - FC ,o- (21) 
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FIG. 4: (Online version in color) Current J as a function 7 for the overdamped case (/i = 0). The 
lines and symbols represent J calculated using the MCFM and MC simulations, respectively. The 
parameters are (a) D x = 1 and Q = 0.2 and (b) D x = 1 and Q = 1, with p = 0.1 (solid lines 
and circles), p = 1 (dotted lines and squares), p = 10 (dashed lines and triangles), and p = 00 
(dot-dashed lines and crosses). 

The current given by Equation [21] does not vanish even when there is no load (i.e. F = 0) 
due to the broken detailed balance. 

In practical calculations of Equation [21], we increased the truncation numbers My. and 
M n in Equation [17] until the current J converged. We also performed MC simulations to 
verify reliability of the MCFM, where the velocity is given by v = [x(T) — x(0)]/T (T = 10 5 ). 
The calculation was repeated 100 times and the average velocity was determined. Below, we 
calculate the dependence of the current J on the relaxation rate 7 (Figured]), the effective 
noise intensity Q (Figure [5]), and the squared variation coefficient p (Figure [6]). 

We first show the current J as a function of the relaxation rate 7 [Equation [3]. In 
Figure H] 7 is plotted against the current J for two sets of parameters: Case (a) D x = 1 and 
Q = 0.2 (Figure [4(a)), and Case (b) D x = 1 and Q = 1 (Figure H](b)). The relationship was 
computed for four values of p (the squared variation coefficient): p = 0.1, 1, 10, and 00 (the 
case p = 00 corresponds to a = and D s = Q/D x ). In Case (a), J became a maximum at 
an intermediate 7, while the current vanished for both 7 — > and 7 — > 00. The small current 
at p = 0.1 was expected because a small p corresponds to weak noise fluctuations. Because 
larger p means larger fluctuations (the kurtosis monotonically increases as a function of p, 
as shown in Figure [3]), it is natural that the current increases as p increases. However, J in 
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FIG. 5: (Online version in color) (a) Current J m; 
case (p = 0). The parameters are D x = 1 with 
(dashed lines), and p = oo (dot-dashed lines). 
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and (b) 7 max as a function Q for the overdamped 
= 0.1 (solid lines), p = 1 (dotted lines), p = 10 
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FIG. 6: (Online version in color) J max as a function of (a) squared variation coefficient p and (b) 
kurtosis n for the overdamped case (p = 0). The parameters are D x = 1 with Q = 0.1 (solid 
lines), Q = 0.2 (dotted lines), Q = 0.4 (dashed lines), Q = 0.8 (dot-dashed lines), and Q = 1.6 
(dot-dot-dashed lines) . 

the range 7 > 10 shows a different tendency and it is higher when p = 1 than when p = 10 
or p = 00. This is due to the noise- intensity modeling in Equation specifically, s(t) in 
Equation [7] for small p rarely has negative values, whereas s(t) for large p has negative as 
well as positive values. For large 7 cases, the "effective" relaxation rate might be measured 
by (|s(t)||s(t')|) rather than by (s(t)s(t')), which indicates that the effective relaxation rates 
for larger p cases are larger than their actual values. In Figure HI the results obtained by 
the MCFM (lines) are always in agreement with MC simulations (symbols), which indicates 
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the reliability of the MCFM. The MC simulations tended to converge for smaller Q cases 
because v is more centered on the mean (thin distributions) when the effective noise intensity 



46[ where noise- 



is smaller. We also note the difference between our model and that of Ref. 
intensity is modulated by a deterministic sinusoidal signal. We found that the width of 
peaks tends to be wider in our model, which is a consequence that the power spectrum of 
Ornstein-Uhlenbeck noise have the Lorentzian function whereas that of a sinusoidal function 
is a delta-peaked function. 

We next show the current J as a function of the effective noise intensity Q [Equation ITT] . 
Since the effect of the relaxation rate 7 varies depending on the squared variation coefficient 
p (as shown above), we removed the 7 dependence of the current J by taking the maximum 
in terms of 7: 

Jmax = Jm ax (Dx,Q,p) = max J(D X , 7, Q, p) , (22) 

7 

7max = lm ax (D x ,Q,p) = argmax J(D X , 7, Q, p) . (23) 

7 

Here, J max and 7 max respectively denote the maximum current and 7 when the current is 
a maximum. Figures [5(a) and (b) show plots of Q against J max and 7 max for D x — 1, 
respectively. Again, computations were performed for four values of p. In Figure [5(a), 
Jmax has maxima as a function Q; the magnitude of this maxima is small for p = 0.1, as 
in Figure [H Transport effects are robust in terms of the noise intensity Q, because the 
current for Q = 2 still exhibits half the maximum current (around Q = 0.2) for large p 
cases. Figure [5(b) shows 7 max (i.e., the relaxation rate yielding the maximum current) as 
a function of Q. For all four values of p, 7 max has a maximum value around Q = 0.1 and 
decreases with increasing Q. This result shows that the noise intensity has to fluctuate over 
a long time scale to enhance the transport capability in noisy environments. 

We also investigated the dependence of p on J max while keeping the other parameters 
constant. Figure [6(a) shows J max (Equation [22]) as a function of p with D x — 1 and five 
values of Q (the effective noise intensity): Q = 0.1,0.2,0.4,0.8, and 1.6. J max increases 
monotonically; the increase exhibits sigmoid-like behavior in terms of log p, indicating that 
■/max is an extremely nonlinear function in terms of p. Figure [6(b) shows the dependence of 
Jmax on the kurtosis k (Equation [13]). Although J max still exhibits nonlinearity as a function 
of n, its nonlinearity is much smaller than that of p. This result shows that the kurtosis can 
be used as an index parameter for Brownian transport driven by noise-intensity fluctuations. 
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FIG. 7: (Online version in color) Current J as a function of p using MC simulations for (a) 7 = 1 
and (b) 7 = 30. The other parameters are D x = 1, Q = 0.2, and F = with p = 1 (circles) and 
p = 00 (squares) . Lines are included as a guide to the eye only. 

3.2. Underdamped case 



We next investigated Brownian transport (Equations and UJ) in the underdamped 
regime, especially from the viewpoint of a mass separation effect [25|. Although we could 
use the MCFM for the overdamped case, the MCFM for an underdamped FPE did not yield 
stable solutions in terms of the current (data not shown) so that we used MC simulations 
only. The underdamped model includes an extra parameter /1 (scaled mass) in addition 
to the four parameters (7, p, Q, and D x ). The velocity is given by v = [x(T) — x(0))/T 
(T = 10 5 ). The calculation was repeated 100 times to obtain the current as the average 
velocity. 

We first show the current J as a function of the scaled mass fi. Figure [7J shows the 
dependence of the current Jon/i for two values of 7 [Case (a) 7 = 1 (Figure [TJa)) and Case 
(b) 7 = 30 (Figure UJb))], where the other parameters are D x = 1, Q = 0.2, F — 0, and 
p = 1 (circles) or p = 00 (squares). In Case (a), the current decreases monotonically and 
the current for p = 00 always exceeds that for p = 1. In contrast, in Case (b), the current 
of p = 00 is always smaller than that of p = 1 and the relation between the magnitudes 
for p = 1 and p = 00 differs from that for Case (a). As shown in Figure |5l the "effective" 
relaxation rate for larger p is larger than the actual relaxation rate. This is the reason why 
converse magnitude relaxation occurs between p = 1 and p = 00 in Case (b). 

We next show the dependence of the average current J on 7 for two values of Q (the 
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FIG. 8: (Online version in color) (a) and (b) Current J with load F > as a function of 7 using 
MC simulations for (a) D x = 1, Q = 0.1, p = 00, and F = 0.025; and (b) D x = 1, Q = 1, p = 00 
and -F = 0.006. J is calculated with MC simulations with p = 0.01 (circles) and p, = 1 (squares), 
and the overdamped case (p = 0) is also shown (solid line), (c) A J (difference of the currents for 
p = 0.01 and 1 defined by Equation I24j) as a function 7 for two mass separations. Lines represent 
LOWESS (locally weighted scatterplot smoothing) of MC simulations (see the text). 

effective noise intensity): Case (a) Q = 0.1 [D x = 1, p = 00, and F = 0.025] (Figure 0(a)) 
and Case (b) Q = 1 [D x = 1, p = 00, and F = 0.006] (Figure 0(b)). In both cases, the load 
F [Case (a) F = 0.025 and Case (b) F = 0.006] increased the ratchet potential for larger 
values of x and the value of F was determined such that particles of p, = 0.01 and p, — 1 
move in opposite directions. Figure shows J of p = 0.01 (circles) and p = 1 (squares) for 
Q = 0.1 (Figure 0(a)) and Q — 1 (Figure 0(b)), where the lines show the MCFM results for 
the overdamped case (p — 0). We include the /i = 0.01 case along with the MCFM result to 
show that results of the underdamped case decrease asymptotically to the overdamped case 
as p — > 0. In Cases (a) and (b), the current J for p = 0.01 has positive values, unlike J for 
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(1=1, which is negative over the whole range of 7. J for (i = 0.01 has peaks at intermediate 
7 values; however, the 7 max values at which the current is a maximum are different in Cases 
(a) and (b): 7 max in Case (a) (Q = 0.1) is around 10, whereas 7 max in Case (b) (Q = 1) is 
about 1. Since the currents for (i = 0.01 and (i = 1 move in the opposite directions, particles 
with (i = 0.01 and (i = 1 can be set apart due to the mass separation effect. To quantify 
the mass separation capability, we define A J as 

AJ = J((i = 0.01) - J((i = 1). (24) 

AJ is a quantity of interest for the separation phenomenon, and systems with larger AJ 
exhibit a better separation capability. FigureEtc) shows A J for two cases [Case (a) (Q = 0.1) 
and Case (b) (Q = 1)] and the lines show locally weighted scatterplot smoothing (LOWESS) 
for each data. It clearly shows that 7 at which the mass separation capability is a maximum 
is different for the two Q cases and that better mass separation is realized at smaller 7 values 
in the larger Q case. This result shows that the mass separation capability is greater for 
slower environmental fluctuations when Brownian motors are subject to strong noise. 



4. DISCUSSION 



We investigated the statistical properties of SIN on Brownian ratchet and found that 
the current was enhanced for smaller 7 in both the overdamped and underdamped regimes 
(Figures [5] and [8]) when the effective noise-intensity Q is larger. This result is intriguing 
because a smaller 7 corresponds to extrinsic fluctuations with a larger time scale. In single- 
cell experiments on E. coli, assuming that mRNA production is Poissonian and that the 
protein burst size has an exponential distribution, the protein copy number obeys the gamma 
distribution 1153: 



™a— 1 „— x/b 

PW = W- <25) 

where a and b are parameters. Observations [15J revealed that the fluctuations (extrinsic 
noise) in a and b are slow and the stationary distribution can be approximated as 

db J o da P(a)P(b), (26) 

where P(a) and P{b) denote the distributions of a and b, respectively. Note that Equa- 
tion [26] is equivalent to the description of superstatistics (see Equation [p. In stochastic 
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gene expression, the intrinsic noise is rather well-explained: its source is stochastic chemical 
reactions with small number of molecules. In contrast, contribution of the extrinsic noise 

IgO ]. Stochastic gene expression is often mod- 



has not been well reasoned or modeled 
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eled with Langevin equations, i.e., continuous approximation of the continuous-time Markov 



chain through Van Kampen's expansion [61] . Many models assume that extrinsic noise affect 
the transcriptional, translational and degradation kinetics as in Equation [261 which results 
in both drift terms and noise-intensity fluctuations of corresponding Langevin equations. 
The stochastic term of Langevin equations is often approximated with the additive white 
noise, and therefore, the extrinsic noise-induced fluctuation fits the overdamped case given 
by Equation and [7J assuming that the drift term fluctuation is negligible compared to that 
of the noise-intensity. Our calculations showed that environmental fluctuations should be 
slow to enhance the transport capability in noisy environments. That is, fluctuations in a 
and h in Equation [25] occur over long time scales. 

Lastly, let us show the correspondence of our calculations with actual time-scales in 
biological experiments. Refs. (54) observed that the time correlation of extrinsic fluctuations 
in E. coli is in the order of cell cycle length T cc . The deterministic part of the protein 
concentration x generally obeys 

where /+ is a protein synthesis term (via translation) and is the degradation rate. In 
Langevin equations, the degradation rate corresponds to the relaxation rate (reciprocal 



of relaxation time) and hence it determines the time scale. Dunlop et al. [55|| defined j3 = 
log2/T cc and concluded that the time correlation of the extrinsic fluctuation is r = T cc / log 2 
(T cc = 60min). Because the ratchet potential given by Equation [S] can be approximated by 
a quadratic function around the minima, Equation [6] can be cast in the form of Equation [271 
with /3 ~ 10. Matching the time scales of these two systems by comparing the relaxation 
time (i.e., comparing (3 in the two cases), the time correlation of the extrinsic fluctuations 
observed in E. coli corresponds to r = 0.1 in our ratchet potential model. Figure [5] shows 
that 7 max is in the approximate range of 1-10, which implies that r max = 7~^ x = 0.1 ~ 1. 
Therefore, the extrinsic noise in the experimental observation is close to optimal for the 
transport effect. Extrinsic noise has non-negligible time-correlation, leaving the system at 
nonequilibrium states. Although there remains a biochemical gap between gene expressions 
and ion transports, our result shows the biochemical advantage of exploiting extrinsic noise 
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for gene regulation. 

In summary, we investigated the transport properties of Brownian ratchet in both over- 



damped and underdamped regimes [44|, |45(. In the overdamped regime, our calculations by 



the MCFM and MC simulations revealed the existence of a maximum current as a function 
of 7 (the relaxation rate) and Q (the effective noise intensity). The maximum current is 
induced at a lower relaxation rate 7 for higher noise intensities. In the underdamped regime, 
MC calculations also showed a maximum for smaller 7 when systems are subject to noisy 
environments. Consequently, the mass separation capability was also maximized for smaller 
7 in such cases. We continue the investigation of ratchet transport in specific biological 
models in future studies. 
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Appendix A: Kurtosis of SIN 

We calculate kurtosis k of SIN defined by 

,= iifgfewn. (A1) 

{{<t)Ut)} 2 ) 

By using an independence relation between s and £ x and the Gaussian nature of £ x , the 
numerator and denominator in Equation IA1I are given by 

({s(t)Ut)V) = 3 (s(t) 4 > (Ut) 2 ) 2 , (A2) 
({s(t)Ut)} 2 ) 2 = (s(t) 2 ) 2 (Ut) 2 ) 2 ■ (A3) 

Let u(t) = s(t) — a. Since u(t) is a standard Ornstein-Uhlenbeck process with (u(t)) = 0, 
the second and fourth-order moments are (u(t) 2 ) = D s and (u(t) 4: ) = 3-D;?, which yields 

(s(t) 2 ) = (u(t) 2 ) +a 2 = D s + a 2 , (A4) 
(s(t) 4 ) = « 4 + 6« 2 (u{tf) + (u{tf) = « 4 + 6a 2 D s + 3D 2 S . (A5) 
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With Equations IA2HA5} the kurtosis is calculated into 



3 (a 4 + 6a 2 D s + 3D 2 ) 
(D s + a 2 ) 



« = ' W (A6) 



Substituting p = D s /a 2 into Equation IA6| we obtain Equation 



Appendix B: Matrix continued fraction method 

We explain the procedure of the MCFM in the overdamped case. Substituting Equa- 
tion [T7] into Equation [151 we obtain a linear algebraic equation: 

n7 + 2F7iki - 4n 2 D x k 2 \ a 2 + 2D S (n+~ 



Cfc,n - ,i. I -- . o i . 

+C fe _i in 7Tfcl + C fc+ i in 7rfci H 7Tfcl H Tiki 

-4C k , n+2 n 2 D x D s k 2 ^(n + 2)(n + 1) - 4C fc)n _ 2 7r 2 J D a!J D a A; 2 V^ - 1) 



-SC^^tHA^VA^ - SC^+^AMVAfa + 1), (Bl) 

where the dimension of the linear algebraic equation ( TBlj) is (2Mfc + l)(M n +l). We show the 
procedure of the MCFM, since a naive conversion of Equation IB 1 1 yields an equation with 
second-nearest-neighbor coupling, which is not compatible with the MCFM (the MCFM 
can treat only first-nearest-neighbor coupling). By introducing = (C^o, C^i ■ • • , Ck,M n ) T , 
Equation IBll is calculated into 



= nklE ^ k ~ 2 + TrifeiEcfc-! + [A - 47r 2 A ! fc 2 B + 2Fvr£;iE] c fc + vrA;iEc fe+1 + 7rfcl ^ Cfc + 2 , ( B 2) 

Here, E is the identity matrix and A and B are (M n + 1) x (M n + 1) matrices defined by 

A„+i,„'+i = -^7<W' ( B3 ) 
B n+ljn /+i = jet 2 + 2D S (^n + J + 2a\/ D a n5 n - 1<n > + 2ayjD 8 (n + l)5„ + i,„', 

+D S v/n(n - l)5 n _ 2 , n / + D s + l)(ra + 2)5 n+2 ,„' , (B4) 



where < n < M n and < n' < M n . Introducing dk = (cj & , cJ fc+1 ) T [58[, Equation lB2l which 
is an equation with second-nearest-neighbor coupling, reduces to the following equation with 
first-nearest-neighbor coupling: 

= Q7/ c fe _i + Q fc c fc + Qfc (B5) 
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where Qfc are (2M n + 2) x (2M n + 2) matrices consisting of submatrices A and B: 

irkiE 2iikiE 
n(2k + l)iE/2 

A - lQn 2 D x k 2 B + AFnkiE 2nkiE 

n(2k + l)iE A - Aw 2 D x (2k + 1) 2 B + 2Frt(2k + l)iE 

7rfciE 
vr(2A; + l)iE 7r(2fc + l)iE/2 
We solve the recurrence relation for Equation IB5I by introducing Sfc and Rfc that satisfy 
Cfc+i = SfcCfc (k > 0) and c^-i = Rfc-iCfc (A; < 0). With Sfc and Rfc, Equation IB5I is calculated 
into 

Sji,-! = - (Q fe + Q^T 1 (fc>l), (B9) 

R fc = - (Q* Rfc-i + Qfc) -1 Qfc (* < -i), (Bio) 

where Sfc and Rfc can be obtained by truncating at large k, namely at k — — Mfc, M&/2). 
For k = 0, we have 

= [Qq R-i + Qo + Q^So] c , (Bll) 

where the first row of the left part of matrix in Equation IB 1 1 1 vanishes due to Equations IB6f 
IB81 Therefore, Co has a nontrivial solution and c"fc can be calculated by recursively applying 
Sfc and Rfc to Cq. 



[1] R. Benzi, A. Sutera, and A. Vulpiani. The mechanism of stochastic resonance. J. Phys. A, 
14:L453, 1981. 

[2] B. McNamara and K. Wiesenfeld. Theory of stochastic resonance. Phys. Rev. A, 39:4854-4869, 
1989. 

[3] P. Jung and P. Hanggi. Amplification of small signals via stochastic resonance. Phys. Rev. A, 
44:8032-8042, 1991. 

[4] L. Gammaitoni, P. Hanggi, P. Jung, and F. Marchesoni. Stochastic resonance. Rev. Mod. 

Phys., 70:223-287, 1998. 
[5] M. D. McDonnell, N. G. Stocks, C. E. M. Pearce, and D. Abbott. Stochastic resonance. 

Cambridge University Press, 2008. 

18 



[6] M. D. McDonnell and D. Abbott. What is stochastic resonance? definitions, misconceptions, 

debates, and its relevance to biology. PLoS Comput. Biol, 5:el000348, 2009. 
[7] F. Marchesoni, L. Gammaitoni, and A. R. Bulsara. Spatiotemporal stochastic resonance in a 

<fi 4 model of kink-antikink nucleation. Phys. Rev. Lett., 76:2609-2612, 1996. 
[8] J. Teramae and D. Tanaka. Robustness of the noise-induced phase synchronization in a general 

class of limit cycle oscillators. Phys. Rev. Lett, 93:204103, 2004. 
[9] J. A. Acebron, L. L. Bonilla, C. J. Perez Vicente, F. Ritort, and R. Spigler. The Kuramoto 

model: A simple paradigm for synchronization phenomena. Rev. Mod. Phys., 77:137-185, Apr 

2005. 

[10] H. Nakao, K. Arai, and Y. Kawamura. Noise-induced synchronization and clustering in en- 
sembles of uncoupled limit-cycle oscillators. Phys. Rev. Lett., 98:184101, 2007. 

[11] M. Kcern, T. C. Elston, W. J. Blake, and J. J. Collins. Stochasticity in gene expression: from 
theories to phenotypes. Nat. Rev., 6:451-464, 2005. 

[12] P. R. Patnaik. External, extrinsic and intrinsic noise in cellular systems: analogies and impli- 
cations for protein synthesis. Biotechnol. Mol. Biol. Rev., 1:121-127, 2006. 

[13] N. Maheshri and E. K. O'Shea. Living with noisy genes: How cells function reliably with 
inherent variability in gene expression. Annu. Rev. Biophys. Biomol. Struct., 36:413-434, 
2007. 

[14] J. Rausenberger, C. Fleck, J. Timmer, and M. Kollmann. Signatures of gene expression noise 

in cellular systems. Prog. Biophys. Mol. Biol., 100:57-66, 2009. 
[15] Y. Taniguchi, P. J. Choi, G.-W. Li, H. Chen, M. Babu, J. Hearn, A. Emili, and X. S. Xie. 

Quantifying E. coli proteome and transcriptome with single-molecule sensitivity in single cells. 

Science, 329:533, 2010. 

[16] J. R. Chabot, J. M. Pedraza, P. Luitel, and A. van Oudenaarden. Stochastic gene expression 
out-of-steady-state in the cyanobacterial circadian clock. Nature, 450:1249-1252, 2007. 

[17] D. Nozaki, D. J. Mar, P. Grigg, and J. J. Collins. Effects of colored noise on stochastic 
resonance in sensory neurons. Phys. Rev. Lett., 82:2402-2405, Mar 1999. 

[18] M. A. Fuentes, R. Toral, and H. S. Wio. Enhancement of stochastic resonance: the role of 
non Gaussian noises. Physica A, 295:114-122, 2001. 

[19] M. O. Magnasco. Forced thermal ratchets. Phys. Rev. Lett, 71:1477-1481, Sep 1993. 

[20] P. Hiinggi and F. Marchesoni. Artificial Brownian motors: Controlling transport on the 

19 



nanoscale. Rev. Mod. Phys., 81:387-442, 2009. 
[21] R. D. Astumian and M. Bier. Fluctuation driven ratchets: molecular motors. Phys. Rev. 

Lett, 72:1766-1769, 1994. 
[22] R. D. Astumian. Thermodynamics and kinetics of a Brownian motor. Science, 276:917-922, 

1997. 

[23] R. D. Astumian and I. Derenyi. Fluctuation driven transport and models of molecular motors 

and pumps. Eur. Biophys. J, 27:474-489, 1998. 
[24] T. Y. Tsong. Na, K-ATPase as a Brownian motor: Electirc field-induced conformational 

fluctuation leads to uphill pumping of cation in the absence of ATP. J. Biol. Phys., 28:309- 

325, 2002. 

[25] F. Marchesoni. Conceptual design of a molecular shuttle. Phys. Lett. A, 237:126-130, 1998. 
[26] H. Linke, T. E. Humphrey, P. E. Lindelof, A. Lofgren, R. Newbury, P. Omling, A. O. Sushkov, 

R. P. Taylor, and H. Xu. Quantum ratchets and quantum heat pumps. Appl. Phys. A, 75:237- 

246, 2002. 

[27] E. Lundh and M. Wallin. Ratchet effect for cold atoms in an optical lattice. Phys. Rev. Lett., 
94:110603, 2005. 

[28] J. L. Mateos. A random walker on a ratchet. Physica A, 351:79-87, 2005. 

[29] S. Bouzat and H. S. Wio. Current and efficiency enhancement in Brownian motors driven by 
non Gaussian noises. Eur. Phys. J. B, 41:97-105, 2004. 

[30] S. E. Mangioni and H. S. Wio. A random walker on a ratchet potential: effect of a non 
Gaussian noise. Eur. Phys. J. B, 61:67-73, 2008. 

[31] E. Prey and K. Kroy. Brownian motion: a paradigm of soft matter and biological physics. 
Ann. Phys., 14:20-50, 2005. 

[32] G. Wilk and Z. Wlodarczyk. Interpretation of the nonextensivity parameter q in some ap- 
plications of Tsallis statistics and Levy distributions. Phys. Rev. Lett., 84:2770-2773, Mar 
2000. 

[33] C. Beck. Dynamical foundations of nonextensive statistical mechanics. Phys. Rev. Lett., 

87:180601, Oct 2001. 
[34] C. Beck and E. G. D. Cohen. Superstatistics. Physica A, 322:267-275, 2003. 
[35] C. Beck. Generalized statistical mechanics for superstatistical systems. Phil. Trans. R. Soc. 

A, 369:453-465, 2011. 

20 



[36] C. Tsallis. Possible generalization of Boltzmann-Gibbs statistics. J. Stat. Phys., 52:479-487, 
1988. 

[37] C. Tsallis. Introduction to Nonextensive Statistical Mechanics: Approaching a Complex World. 
Springer, 2009. 

[38] S. M. D. Queiros and C. Tsallis. On the connection between financial processes with stochastic 
volatility and nonextensive statistical mechanics. Eur. Phys. J. B, 48:139-148, 2005. 

[39] C. Beck. Superstatistical Brownian motion. Prog. Theor. Phys. Suppl, 162:29-36, 2006. 

[40] R. F. Rodriguez and I. Santamaria-Holek. Superstatistics of Brownian motion: A comparative 
study. Physica A, 385:456-464, 2007. 

[41] P. Jizba and H. Kleinert. Superpositions of probability distributions. Phys. Rev. E, 78:031122, 
2008. 

[42] Y. Hasegawa and M. Arita. Bistable stochastic processes in the g-exponential family. Physica 
A, 389:4450-4461, 2010. 

[43] G. C. Yalcin and C. Beck. Currents in complex polymers: An example of superstatistics for 

short time series. Phys. Lett. A, 376:2344-2347, 2012. 
[44] Y. Hasegawa and M. Arita. Noise-intensity fluctuation in Langevin model and its higher-order 

Fokker-Planck equation. Physica A, 390:1051-1063, 2011. 
[45] Y. Hasegawa and M. Arita. Escape process and stochastic resonance under noise-intensity 

fluctuation. Phys. Lett. A, 375:3450-3458, 2011. 
[46] P. Reimann, R. Bartussek, R. Haufiler, and P. Hanggi. Brownian motors driven by temperature 

oscillations. Phys. Lett. A, 215:26-31, 1996. 
[47] Y.-X. Li. Transport generated by fluctuating temperature. Physica A, 238:245-251, 1997. 
[48] Y. Zhang and J. Chen. Investigation on a temporal asymmetric oscillating temperature 

ratchet. Physica A, 387:3443-3448, 2008. 
[49] M. Borromeo, S. Giusepponi, and F. Marchesoni. Recycled noise rectification: An automated 

Maxwell's deamon. Phys. Rev. E, 74:031121, 2006. 
[50] W. A. M. Morgado, S. M. D. Queiros, and D. O. Soares-Pinto. On exact time averages of a 

massive Poisson particle. J. Stat. Mech., page P06010, 2011. 
[51] C. R. Doering and J. C. Gadoua. Resonant activation over a fluctuating barrier. Phys. Rev. 

Lett, 69:2318-2321, Oct 1992. 
[52] R. N. Mantegna and B. Spagnolo. Noise enhanced stability in an unstable system. Phys. Rev. 

21 



Lett., 76:563-566, Jan 1996. 
[53] B. Spagnolo, N. V. Agudov, and A. A. Dubkov. Noise enhanced stability. Acta Phys. Pol. B, 
35:1419-1436, 2004. 

[54] N. Rosenfeld, J. W. Young, U. Alon, P. S. Swain, and M. B. Elowitz. Gene regulation at the 

single-cell level. Science, 307:1962-1965, 2005. 
[55] M. J. Dunlop, R. S. C. Ill, J. H. Levine, R. M. Murray, and M. B. Elowitz. Regulatory activity 

revealed by dynamic correlations in gene expression noise. Nature Genetics, 40:14393-1498, 

December 2008. 

[56] R. Bartussek, P. Reimann, and P. Hanggi. Precise numerics versus theory for correlation 

ratchets. Phys. Rev. Lett, 76:1166-1169, 1996. 
[57] B. Lindner, L. Schimansky-Geier, P. Reimann, P. Hanggi, and M. Nagaoka. Inertia ratchets: 

A numerical study versus theory. Phys. Rev. E, 59:1417-1424, 1999. 
[58] H. Risken. The Fokker-Planck Equation: Methods of Solution and Applications. Springer, 

2nd edition, 1989. 

[59] M. Scott, B. Ingalls, and M. Kaern. Estimations of intrinsic and extrinsic noise in models of 

nonlinear genetic networks. Chaos, 16:026107, 2006. 
[60] V. Shahrezaei, J. F. Ollivier, and P. S. Swain. Colored extrinsic fluctuations and stochastic 

gene expression. Mol. Syst. Biol, 4:196, 2008. 
[61] V. Kampen. Stochastic Process Theory in Physics and Chemistry. North-Holland, 1992. 



22 



